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The energy and zenith angle distribution of neutrino induced, upward going muons 
can give direct information on the presence of i^-oscillations in precisely the range of 
parameters suggested as a solution of the atmospheric neutrino problem. We discuss here 
the uncertainties in the theoretical prediction. The shape of the zenith angle distribution 
of the muon flux is quite insensitive to modifications of the theoretical input and is a 
good probe for the existence of neutrino oscillations. We conclude that the existing data 
sample on (/-induced muons has the statistical power to confirm or refute the //-oscillation 
solution of the atmospheric neutrino problem. 

PACS numbers: 14.60. Pq, 96.40.Tv 



^ I The recent results of the Superkamiokande (SK) experiment 0| on atmospheric neutrinos confirm the 

00 . existence of a problem in the relative numbers of Ve and induced events, that may be interpreted 
as an indication of neutrino oscillations with large mixing. The favoured Am^ has however been shifted 
, to values lower than those suggested by the older Kamiokande results Q. The results of the Chooz 
experiment Q| exclude large amplitude oscillations of electron antineutrinos into other antineutrinos in 
, the region of interest for Ato^, however v^^ <-> Vr (or Vs where Vs is a sterile neutrino) oscillations 

■ remain a viable solution of the problem. A dedicated long-baseline experiment could provide a clear 
confirmation of the existence of oscillations ||^ ; the actual detection of r-leptons produced in the detector 
would be especially striking, however this would only be possible with a high energy neutrino beam and 
the distance needed to observe such an effect would become prohibitively long if Am^ 10~'^ eV^. It is 
Q-i. therefore very desirable to obtain clearer indications on the value of the parameter Am^ . 

• In this work, we rediscuss the information that can be obtained from the analysis of z/-induced upward 

going muons. In particular we consider in detail the shape of the zenith angle distribution and suggest 
that it is a sensitive and model independent probe of the existence of neutrino oscillations. In the presence 
of r/-oscillations the muon flux is reduced in its absolute value, its energy spectrum is distorted, and the 
zenith angle distribution is deformed. Since the v induced flux can be measured only in half of the solid 
' angle, and the deviations from the up-down symmetry (j)^{E,9) = 4>fi{E,'K — 9), valid in the absence of 
oscillations, cannot be measured, the interpretation of the data requires a comparison with a calculated 
prediction, and the estimate of theoretical uncertainties in the calculation is crucial. 

The absolute value of the predicted flux is affected by a large uncertainty ^ 20% This uncertainty 
is not relevant in the study of the shape of the muon energy spectrum. Since no underground muon 
detector has used a magnetic field, a measurement of the energy can only be obtained for the muons that 
stop in the detector using the particle range. In the presence of oscillations the flux of low energy muons 
is suppressed more, and the ratio r — Ns/Np of the numbers of 'stopping' and 'passing' muon events is 
reduced with respect to the no-oscillation prediction. From a measurement of r, the 1MB collaboration 
has obtained an exclusion region in the (Am^, sin^ 2 0^,-) plane that is in sharp contradiction with the 
allowed region obtained by SK [Q from the analysis of the events with v interactions inside the detector. 
The significance of the discrepancy is much reduced in the critical reanalysis of , where the prediction 
for the value of r in the absence of oscillations is increased by ~ 12% with respect to the calculation used 
in the original 1MB work. The SK detector, with a larger effective area, larger containement volume and 
better spatial and energy resolution can provide a new and more accurate measurement of the stop/pass 
ratio. 

The direction of the muons is well measured in underground detectors and the shape of the angular 
distribution can be studied in detail and compared to different predictions (with or without oscillations) 
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using statistical tests. For illustration and the sake of clarity in this work we will simply use the ratio of 
the of 'horizontal' and 'vertical' muon events, a single number, to describe the shape of the distribution: 
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In the presence of oscillations the flux of vertical muons is supressed more (because the quasi-collinear 
parent neutrinos have longer pathlengths) and S increases, by as much as 30% for Am^ in the region 
of interest and maximal mixing. The separate measurement of the quantities Sg and Sp, related to the 
angular distributions for 'stopping' and 'passing' muons, results in an improvement of the sensitivity, a 
better control of systematic effects, and in case of a positive signal a more precise indication of the value 
of Am^. We estimate the uncertainty in the predicted value of S, Sg and Sp as ^ 4%. 

In previous theoretical discussions of stopping muons it was assumed for simplicity that stopping 
muons are in the energy interval -Emin < Ep < E* , and passing muons have energies Ep > E*. This 
approximation is clearly not realistic. For simple geometrical reasons, the energy interval where stopping 
muons can be identified depends on the zenith angle; moreover in general, even for a fixed zenith angle, 
muons of the same energy can either stop in or pass through the detector, depending on their trajectory. 
In general the rates of stopping and passing muons in the detector can be estimated as: 



dcos 9 
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where 4>p{Ep, cos 0) is the differential flux of upward going muons of energy and zenith angle 9, 
and Ag(^p-j(Ep,9) is the detector acceptance for stopping (passing) muons. A detailed calculation of the 
acceptances as a function of muon energy and direction is only possible for the experimental collaboration. 
In the following, considering the case of the SK detector, as a first approximation we describe it as a 
homogeneous cylindrical volume of radius R = 16.9 m and height H = 36.2 m and approximate 
the muon trajectories as straight lines with length equal to the muon range in water; in this way the 
calculation of the acceptance is reduced to an elementary geometry problem. 

The requirement that the muon track be entirely contained in the detector volume gives: 

ME^,9) = A{LiE^),9), (3) 

where L{E) is the muon range in water and 

A{L,9) ^2 RH sin 9 ^/T^^ +2 R^\cos9\ [cos-^ x - 3 x ^/T^~^] e[L„^ax{0) - L] , (4) 

with X — Lsm9/2R and imax(^^) = min[2i?/ sm9, iJ/|cos6'|]. The acceptance for stopping muons is 
non-zero only when the muon range is in the interval imin < E{E) < L*, where Lmin is the minimum 
track length required for detection (7 m in the first analysis of SK and L* is the maximum linear 

dimension of the detector (49.5 m in SK); this corresponds to the energy interval 1.45 < Ep < 10.63 GeV. 
The acceptance for passing muons is given by: 

Ap{Ep,9) ^ A{Ln,in,0) - A,{Ep,9) . (5) 

In Table I we give our estimates of some quantities of interest (total event rate, r, Ss, Sp) calculated 
in the absence of oscillations with different choices of the theoretical input, as described in the first three 
columns. Column 1 indicates the choice of a neutrino flux from those available in the literature [pT|-p7t. 
In column 2 the label LLS indicates that the contributions of quasi-elastic scattering and of single pion 
production to the neutrino cross sections have been explicitely included |Q , while the label DIS indicates 
that the neutrino cross section is described by the formulae of deep inelastic scattering; the more accurate 
treatment of the lowest multiplicity channels in the LLS option results in a higher total flux and a higher 
stop/pass ratio r. The third column indicates the choice made for the Parton Distribution Functions 
(PDF), using the label defined in the library code PDFLIB ||2^]. The first row in table I contains the 
absolute values of the predictions for our reference model (Bartol (pi, , PDF by |Q and LLS § cross 
sections); in the other rows we report ratios with respect to the reference values. 
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Our reference model predicts a rate of 521 yr ^ of passing muons with track length larger than 7 m 
This can be compared with a preliminary result obtained in 229 days of live time of 425 ± 26 yr^^ 10 



We also predict a rate of 218 yr~^ stopping muons. It is important to observe that this prediction is quite 
sensitive to the value of the cut iminj because the muon flux is steeply falling with energy in this region; 
a lower (higher) value Lmin = 5 (9) m results in a stopping rate of 269 (169) yr~^. It will be therefore 
important to consider in detail the resolution in the measurement of the stopping point of the muon 
tracks. The distribution of the length L of the stopping tracks is a measurement of the muon spectrum 
in the low energy region; however this region is narrow, and with the expected event rates the possible 
distortions of the dN/dL distributions will be hard to detect. 

Our predictions for the parameters Ss and Sp that describe the shape of the zenith angle distributions 
of stopping and passing muons are 1.45 and 1.63. In both cases there is an excess of horizontal muons, 
and the anisotropy is stronger for the higher energy (passing) muons, Sp > S^- The anisotropy of the 
muon flux and its energy dependence originate from the zenith angle dependence of the neutrino flux. 
The angular distribution is approximately isotropic for low energy neutrinos, and it develops a stronger 
and stronger dependence on the zenith angle with increasing neutrino energy. This can be understood by 
considering the larger decay probability for high energy mesons and muons traveling in a direction near 
to the horizontal. 

The spread of numbers in the columns of table I can be used as an estimate of the uncertainty in the 
prediction for the different observables. We estimate Ar ~ 8%, ASs ~ 3 % and ASp ^ 4 %. The MMK 
neutrino flux gives 5's(p) larger than our reference value by a factor 1.16 (1.13), in contrast to all 
other models, and we disregard this prediction in the following. 

To explore the possible effects of modifications in the theoretical inputs of the neutrino flux calculation, 
we have repeated the i/-fiux calculation of p^ ] including some extreme changes in the initial assumptions. 
In [p^ the primary cosmic ray flux has a power law energy spectrum E~'^ with a = 2.7; steepening 
(flattening) the spectrum using a — 2.8 (2.6) increases (decreases) the stop/pass ratio r by +26% (- 
21%), correspondingly the cos 6* distribution becomes flatter (steeper) with the much smaller variation 
of =f1% in Ss and +3% in Sp. Increasing (reducing) the K meson yield by a factor 1.5 changes r by 
±8%, Ss by ±3% and Sp by ±4%. These tests indicate that the prediction on the shape of the angular 
distribution is quite stable. 

The effects of neutrino oscillations on the observables that we are discussing are illustrated in the 
figures. In fig. 1 the solid curve shows the dependence on Am^ for maximal mixing of the ratio r 
between the number of stopping and passing events, for our reference model. The error bars indicate 
only statistical errors estimated for two years of SK data-taking (i.e. in a few months from now). In 
the region 10""^ < ATO^(eV^) < 10~^, favoured by the anomaly in the contained and semi-contained 
events, the stop/pass ratio has a value r < 0.30 to be compared with a no-oscillation prediction 0.42. 
It is apparent that in this region of Am^, for large mixing, the statistical significance of the effects of 
oscillations is clear (as large as ^ 7 u), and the theoretical uncertainty of the prediction is likely to be 
the most important source of error. With the estimate of 8% for the theoretical error and reasonable 
systematic errors in the experiment, an effect should be clearly visible if neutrino oscillations are indeed 
the reason of the atmospheric neutrino anomaly. However, the suppression of r is approximately the 
same for Am^ = 10""^ or 10^^ eV^ and therefore with an analysis of this quantity it will not be possible 
to resolve the present uncertainty on Am^ . In fig. 1 we also included (dashed line) the Am^ dependence 
for maximal mixing of r for oscillations into sterile neutrinos, the curve being shifted because of the 
matter effects p. 

We have plotted in fig. 2 versus Arri^, again for maximal mixing, the quantities Ss and Sp that are a 
measure of the shape of the angular distributions. The error bars indicate the estimated statistical errors, 
which in this case are likely to be the dominant source of error. In the presence of neutrino oscillations 
both Ss and Sp are larger than the no-oscillation expectation because of the larger suppression of the 
vertical flux. For a fixed value of the mixing parameter, the effect of oscillations on Ss (Sp) has a broad 
maximum for Am^ ~ 0.7 • 10"'^ (10~^) eV^. Qualitatively we can expect that the angular distortion 
produced by i^-oscillations is maximum when 

^ ^ ^E,) , (6) 
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where i?^ is the earth radius and {E^) is the typical energy of the neutrinos that produce the upward 
going muon signal; for smaller values of Aw? the oscillations do not have time to develop, for larger values 
the effects of oscillations suppress the flux equally for all directions. The signals of passing and stopping 
muons are produced by neutrinos with median energy {Et^^^) ~ 8(100) GeV and this is reflected in the 
positions of the maxima in fig. 2. We also note that for large mixing and Am^ ~ 0.2-1 • 10"'^ eV^ one 
should have Ss > Sp, an essentially model-independent signal of oscillations, albeit with low statistical 
significance. The matter effects, present in the i/^ ^ Vs case, induce a detectable but different deformation 
in the angular distribution; the parameters Ss^p are not a good measure of the effect in this case, since 
the suppression for 'horizontal' and 'vertical' events are approximately equal, after integration over the 
solid angle. 

In fig. 3 we present in the usual (Arn^, sin^ 2 0^1-) plane the regions that may be excluded at 95% c.l. 
after two years of SK running assuming that the experiment measures the central value of the expectations 
(without neutrino oscillations), by a measurement of r (solid fine), Ss (dashed) and Sp (dot-dashed). In 
order to produce this plot, we (arbitrarily) combined quadratically the expected statistical error and 
the theoretical uncertainty. Conversely, the effect of atmospheric i/^ <->• Ur oscillations for representative 
parameter values, maximal mixing and An? = 10"'^ (10~^) eV^, would result in a reduction in the 
stop/pass ratio r of 30 (26)%, with a deviation of 3.2 (2.7) cr's including the theoretical error, and shape 
parameters Ss = 2.14 ± 0.27 (1.55 ± 0.21) and Sp = 1.82 ± 0.12 (2.13 ± 0.17), that correspond to 2.5 
and 1.6 (0.5 and 2.9) standard deviation effects. In both cases the signal should be detectable, and the 
analysis offers a rough measurement of the Am? parameter. 

Experimental results on upward going muons have already been presented. Baksan [T^ , MACRO 
||l9[] , Kamiokande j2^, 1MB |2^,01, and Superkamiokande [0, have collected large statistical samples of 
'passing' events. The results for the total rates are inconclusive, with respect to the oscillations, because 
of the large errors and theoretical uncertainties. All these experiments have presented measurements 
of the muon flux above a threshold energy as a function of the zenith angle. We point out that the 
extraction of a detector independent flux from the measured passing muon rates is both non trivial 
and model dependent, because no detector has a well defined and angle independent threshold, and one 
needs corrections that make use of the theoretical prediction. Since this is normally done assuming the 
absence of oscillations, the fluxes can be used to test this hypothesis, but great care has to be taken 
in the extraction of oscillation parameters. The comparison of event rates with detailed MonteCarlo 
predictions is free from these interpretation problems, that is why we preferred to make a rough estimate 
of the acceptance and discuss predictions of directly measurable quantities. For similar reasons, it seems 
dangerous to reduce the data of different experiments to a common threshold, as was done in a recent 
analysis ||2^. The experimental results are summarized in table II, where we give the ratio S^^^ of the 
vertical and horizontal fluxes obtained from the published data, and our no-oscillation prediction S^. 
All experiments measure a value larger than the prediction. This could be interpreted as an effect of 
neutrino oscillations, but such a conclusion would be premature. The fit of the shape of the angular 
distributions for the individual experiments is poor even after the inclusion of oscillations; moreover a 
comparison between experiments shows that the normalizations and detailed shapes of the data are not 
in good agreement with each other, suggesting the presence of unforeseen systematic effects (see Q for 
a discussion) . Nonetheless the qualitative result of an excess of horizontal muons in all experiments is a 
hint that should be studied in more detail. 

In summary, a better measurement of upward going muons, as it will be possible in the Su- 
perKamiokande detector, can test the neutrino oscillation solution of the (semi)contained atmospheric 
neutrino problem. Sensitive quantities are: the ratio of stopping to passing muons, already noted and 
discussed, and the shapes of angular distributions, that are affected by an even smaller theoretical un- 
certainty. 
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TABLE I. Predictions for observables related to upward going muon flux in the absence of oscillations for the 
SuperKamiokande detector. 
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r 


Ss 


Sp 


Bartol 


LLS 


GRV94LO 


738 


0.42 


1.45 


1.63 


HKKM 


LLS 


GRV94LO 


0.94 


0.98 


1.00 


1.00 


BDZ 
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GRV94LO 
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1.00 
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Volkova 
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FLUKA 
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1.01 


1.02 


MMK 


LLS 


GRV94LO 


1.00 


1.09 


1.16 


1.13 


Bartol 


LLS 


GRVLO 
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1.00 


1.00 


1.00 


Bartol 
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GRV94DIS 


0.98 


0.99 


1.00 


1.00 
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CTEQ4-23 


1.05 


0.98 


1.00 


1.00 


Bartol 


LLS 


EHLQ2 


0.90 


1.02 


1.00 


1.00 


Bartol 


DIS 


GRV94LO 


0.99 


0.96 


1.00 


1.00 


Bartol 


DIS 


GRVLO 


0.92 


0.95 


1.00 


1.00 


Bartol 


DIS 


EHLQ2 


0.85 


0.89 


1.00 


1.00 



TABLE II. Measurements of the shape of the zenith an- 
gle distribution of the i^-induced muon flux. The theoretical 
expectation is based on our reference model. 



Detector 


Et (GeV) 


riexp 


nth 


Baksan 


1.0 


1.82 ±0.18 


1.50 


MACRO* 


1.0 


2.97 ±0.65 


1.50 


1MB 


1.8 


1.79 ±0.18 


1.53 


Kamiokande 


3.0 


1.91 ±0.21 


1.56 


SK 


6.0 


1.83 ± 0.23 


1.60 



^The Macro detector has a very low acceptance for near to horizontal 
muons and the datum in the first angular bin has a very large error, 
that propagates into the ratio of fluxes. If one neglects the first bin 
and extrapolates from the others a ratio similar to the other rows 
is obtained. 
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FIG. 1. Ratio of the rates of 'stopping' and 'passing' muons (with Lmin = 7 m) in SuperKamiokande for maximal mixing. 
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FIG. 2. Ratio of the 'horizontal' and 'vertical' muon rates for stopping and passing events in SuperKamiokande for v^-Vr 
oscillations with maximal mixing 
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FIG. 3. Regions (to the right of the corresponding lines) in the (Am^,sin^ 2S^t) plane, that can be excluded at 95% C.L. 

with two years of SupcrKamiokandc data from incasurcrncrits of r, Ss ans Sp. The exclusion curves are calculated assuming 
absence of oscillations and including theoretical errors of 8, 3 and 4% for the predictions. 
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